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Abstract 

We study the semileptonic decays of heavy mesons into hght pseudoscalars by 
making use of dispersion relations. Constraints from heavy quark symmetry, chiral 
symmetry and perturbative QCD are implemented into a dispersive model for the 
form-factors. Large deviations from i?*-pole dominance are observed in S — > niv. We 
discuss the model prediction for this mode and its possible impact on the extraction 

of iKfel- 
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1 Introduction 



Semileptonic decays of heavy hadrons are of great interest given that they are the cleanest 
way to probe the mixing between quarks of the third generation with those of the hghter 
famihes. The extraction of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element Vet 
from the b — > ciu transitions is largely freed from theoretical uncertainties since the advent 
of the heavy quark symmetries This is the case for exclusive decays as well as for 
the inclusive lepton spectrum . On the other hand, the b —>■ uiu transitions involving the 
CKM element Vub are still plagued with large theoretical uncertainties. The inclusive lepton 
spectrum above the b ciu end-point, from which can be extracted, is still the 

realm of models due to the breakdown of the heavy quark operator product expansion. For 
the exclusive decays, the use of heavy quark symmetry (HQS) is reduced to relating the D 
and B decays to light hadrons at fixed values of the recoil energy. The main shortcoming 
of this prediction for B decays is that it only covers recoil energies available in D decays. 
This is particularly troublesome in the vr mode, for which most of the rate might be at large 
recoil energies. It would thus be very useful to have a full calculation of the B Tciu decay 
consistent with the constraints from HQS as well as chiral symmetry for heavy hadrons. 
The information from HQS is contained in the scaling behavior of the form-factors with 
the heavy masses as well as in the spin symmetry relations for members of the same HQS 
spin multiplet. In the decays H Triu, Chiral Perturbation Theory for Heavy Hadrons 
(CPTHH) IQ, 1^ promotes nearest singularity (pole) dominance to the leading contribution 
in the chiral and heavy quark expansions when the pion is soft. The validity of nearest pole 
dominance in the soft pion limit is not a surprising result given the proximity of the pole 
to the physical region where the pion recoil is small. However, there is no reason to believe 
that this is also true at higher recoil momenta. Deviations from the pure pole behavior at 
pion recoil energies above ~ IGeV in B Triu result in large modifications in the branching 
ratio. Conversely, the physical region in the case of the D modes is confined to be closer to 
the pole and therefore one expects the approximation of the form-factors by the single pole 
to be a rather good one over a large fraction of phase space. In this paper, we construct a 
model which reflects pole dominance at low pion energies but at the same time includes other 
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effects that are potentially important, mostly in B decays. These will include the effect of 
resonances other than the i?*-pole as well as the continuum. A natural model emerges from 



determined not only by the isolated poles but also by the singularities above the {uih + m.„y 
threshold. The contributions from resonances above threshold can be parametrized in a way 
compatible with HQS. This allows us to fix the parameters of the model in D decays and 
have a prediction for all the available phase space in B —>■ -kIv. 

Recently CLEO has observed this decay for the first time and measured its branching 
ratio [[^]. It is therefore imperative to realistically assess the theoretical uncertainties asso- 
ciated with this mode and their impact on the extraction of \Vub\- The model we present 
here is an attempt to address this issue. We expect the contributions to the dispersion re- 
lations governing these transitions to be highly constrained by chiral symmetry, HQS and 
perturbative QCD, to the point of having a complete picture of the physics involved and very 
little freedom. In the next section we discuss the various contributions to the form-factors in 
H — > vr^z/, with H = {D, B), in the language of dispersion relations. In Section 3 we derive 
constraints from chiral symmetry, HQS and perturbative QCD which will highly determine 
those contributions. In Section 4 we present a model that naturally emerges from all the 
theoretical constraints, discuss its predictions and compare them with other calculations. 
We summarize our results and conclude in Section 5. 

2 Dispersion Relations 

The hadronic matrix element for the B^ TT^i^u transition can be written as 



where = {P — p.,^)^ is the momentum transferred to the leptons. In the approximation 
where the leptons are massless, only the form-factor f+{q^) enters the partial rate. This 
form-factor obeys a dispersion relation of the form 



the dispersion relations for the form-factors 0. Their properties in the physical region are 



(7r(pJ|M7^6|/J(P)) = U{q^){P + p^), + f.{q^){P-p^) 



(1) 




(2) 
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where sq = {niH + ^vr)^- The isolated pole at mn* < sq is actually present in the B meson 
case, whereas for the D ^ it transition the D* pole is located almost exactly at threshold. 
There are no anomalous thresholds in the H —>■ iriu transitions, although there might be 
in H ^ piu. Contributing to the imaginary part in (H) are all possible intermediate states 
that couple to Hn and are annihilated by the weak vertex, as shown schematically in Fig. 1. 
These include the multiparticle continuum as well as resonances. The latter must be radially 
excited = 1^ states in order to contribute to f+(t), and are expected to be located 
at about ^ 1 GeV above the ground state [§. The first contribution appearing above Sq 
corresponds to the H^*^tx continuum. Other contributions to the continuum include states 
with H'^*'> mesons and various light mesons (e.g. multipion states). We will neglect them at 
these values of s because they involve higher order terms in the chiral expansion. Finally, 
we consider the contributions of orbitally excited H mesons with one pion. In principle 
these could be important given that the L = 1, P-wave states are only about 500 MeV 
above the ground state [^]. The lightest L = 1 states correspond to (0^, 1^) and (l+,2+) 
doublets. The second doublet, however, does not couple to the ground state to leading order 
in CPTHH The doublet (0^, 1^) does couple to the ground state doublet, but the vertex 



Hn — > (0^, l+)7r vanishes to leading order. Therefore, not far from threshold, the continuum 
contribution to the dispersive integral can be approximated by the H^*^7c intermediate 
states. We will estimate this contribution in CPTHH in Section 3.1. 

The contributions from the radial excitations of the H* dominate the imaginary part at 
intermediate values of s. Separating them from the continuum one can express (0) as 

hit) = ,\, +- / \ .[ +j:<^^Mt). (3) 

{mjj* -t) IT J so [s - t - Id) 

where the cutoff in the continuum integral defines the maximum center of mass energy in Hit 
scattering for which the main contribution comes from the continuum. This also corresponds 
to the beginning of the resonance-dominated region. Typically, A defines an energy about 
0.7 to 1 GeV above threshold. Therefore, the use of CPTHH to compute the continuum 
contribution to the dispersive integral is justified. 
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In the narrow width approximation 



1 



(4) 



t 



whereas if finite width effects are taken into account these functions take the form 



7li{t) = — i arctan 



(5) 



with Mi the mass of the i-th radial excitation. In deriving the widths Fj were assumed 
to be constant. The a^'s are the couplings analogous to 7 of the nearest resonance if*, the 
residues at the poles. They involve the strong coupling between the z-th resonance and if vr, 
as well as the decay constant of the resonance. 

The physical region for the decay H — * vrfi/ is given by the interval t = (0,tniax), with 
^max = i^H — ^n)'^- Thus the if*-pole, the first term in corresponds to the singularity 
closest to this region and will be the dominant contribution for values of t close to tmax- 
However, the question of how good this approximation is in each case is not a simple one. For 
instance, neglecting the continuum contribution, it is likely that the if *-pole approximation 
will be a good one as long as the three- momentum of the recoiling pion, p^, is not larger 
than Ai = Mi — mjj{,), the gap between the ground state and the first radially excited 
state. This intuitive picture suggests that the if *-pole term in should be a reasonable 
approximation to for D semileptonic decays, given that almost all its phase-space falls 
in this region. However, this is certainly not the case for B iriu. Although the spacings 
between resonances and the ground state are independent of the heavy quark mass to a very 
good approximation, now the recoil momentum of the pion can be as large as p™^^ ^7713/2. 
As increases and we move away from the if *-pole, the relative infiuence of the higher 
resonances grows. These deviations from the pure if *-pole behavior at large values of p^^ are 
particularly important given that the pion momentum distribution goes as p^. As we will 
see below, large changes in the total rate and the shape of the t distribution occur when the 
corrections to the if *-pole behavior at large p.^ are taken into account. 

In what follows, we will consider the theoretical constraints that can be imposed on f+{t). 
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These constitute the basis for a model calculation of the B ^ it transition form-factor, which 



will incorporate all of these constraints. 

3 Theoretical Constraints 

Although the form-factor /+(t) is not calculable in perturbation theory, it must satisfy several 
constraints. These result from the application of HQS, chiral symmetry and the asymptotic 
behavior imposed by perturbative QCD for exclusive processes. 

3.1 Chiral Symmetry 

Chiral Perturbation Theory for Heavy Hadrons (CPTHH) provides a formal framework for 
the approximation that only keeps the first term in (|^) @, 3. In the effective theory that 
couples heavy hadrons to goldstone bosons respecting HQS and chiral symmetry, the heavy 
meson fields are represented by the 4x4 matrices 

n = - Hi,} (6) 

where and H are the 1~ and 0~ ground state fileds respectively, and is the heavy 
meson four-velocity. The goldstone bosons enter through 

^ = exp{iTT{x)aTa/ f) (7) 

with Ta, (a = 1, . . . , A^^ — 1) the SU{N) generators and na(x) the goldstone boson fields. 
To leading order f = f-jr, the pion decay constant. The leading order Lagrangian, invariant 
under HQS and chiral symmetry, is given by @, |^ 



£efr. = iTr 



Hv.DH 



+ gTr n^ni5 



The coupling constant g is independent of the heavy mass and the axial- vector field is defined 
by 
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Requiring that the weak current transforms as a left-handed doublet implies 



(10) 



with fn the heavy pseudoscalar decay constant. This completes the description to leading 
order both in l/rriH and l/A^^,, where is the scale of chiral symmetry breaking. The 
H — s> Tciu transition receives a direct contribution from ([l^), as well as an H*-pole term 
resulting from the H*H7t interaction governed by the coupling g in (j^) followed by the H*- 
vacuum transition governed by (|10|). The latter dominates the former in the l/rriH expansion. 
Thus, for soft pion momentum, the transition form-factor can be approximated by 

Therefore, CPTHH tells us that to leading order the form-factor is approximated by the first 
term in the dispersion relation (|^) with 

7= ^ ■ (12) 

In this way CPTHH provides an approximate normalization of /+(t) at low pion recoil 
momentum, ^ IGeV. 

On the other hand, CPTHH can be used to compute the continuum contribution in (|^) 
for values of the integrand close to threshold. At this values of s unitarity implies 

lm\rr'\s)\ = a{s)T\s)j^{s) (13) 

where the threshold factor is 

a{s) ^ (l - ^^^L^^Iht^ [x - K,±Z!^ j (14) 

and T(s) is the i/vr Hti scattering amplitude projected onto the J = 1 partial wave. We 
neglect the contribution of the if*7r intermediate state which is suppressed by a factor of 
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P-k/^h relative to ifvr. Therefore, above sq but below the resonance region in the cut, the 
phase of f+{s) is given by the -ffvr Htc scattering phase-shift 

sin5+(s) = |T(s)|. (15) 

The computation from (|^) is straightforward. After projecting to the correct partial wave 
and isospin channels the phase, to leading order, is given by 

2Att\UJ \v.p^-A v.p^ + AJ 
where A = mn* — ttlh is the mass splitting within the ground state and 

s — mir — ml 

^•P- = 1^ 

z mil 

is the pion energy in the H rest frame. The dominant effect, once again, comes from the 
-?/*-pole exchange, both in the s and t channels. These are governed by the same coupling 
g entering in the H*Ht{ vertex in (|^), and therefore do not introduce any new parameter. 
The phase in (|T^) together with the leading order expression for |/+(s)| yields /m[/™'^*-](s), 
in the threshold region. The corresponding dispersive part provides, in a next-to-leading 
singularity approximation, the first correction to the if*-pole behavior. However, as we will 
see below, this is not the most important modification to pole dominance coming from the 
cut. The presence of radially excited states, coupling to both the ground state and the weak 
current turns out to be a more significant correction. 

3.2 Heavy Quark Symmetry 

Although HQS is an ingredient of the CPTHH formulation, the applicability of the flavor 
symmetry is very limited in practice, as it was mentioned above, if the scaling of f+{t) with 
mil is used to relate D -niv to B ^ Tiiu. However, we will show here that the application 
of HQS to the dispersion relation of (0) does not suffer from such limitations. The first two 
terms in (0) already have a defined behavior with mn built in CPTHH. On the other hand. 
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the residues a, governing the resonant contributions scale as m^:^ , the same way 7 does. 
This allows us to write 



where gi is a dimensionless constant characterizing the coupling of the radially excited reso- 
nance Hi to H^*^TT, and is essentially the ratio of the mass independent decay constant of 
Hi to that of the ground state, fn- On the other hand, the spacing among resonances and 
also the gap to the ground state is independent of ttih, to leading order in the heavy mass. 
Therefore, if we knew the a^'s and the masses of the resonances for the charm system, we 
would know for B —>■ Triu in the whole physical region. This is an important departure 
from the application of the flavor HQS to semileptonic decays. The HQS is applied to the 
resonances in the cut, which have excitation energies independent of the heavy mass as long 
as both the c and the b quark are considered to be sufficiently heavy. Of course, we do not 
know a priori the values of the a^'s and the spacings among resonances and the ground state. 
The latter can be calculated in potential models that have been successful in predicting the 
spectrum of the orbitally excited heavy mesons [^j and are expected to yield a good approx- 
imation also for the radial excitations. Regarding the couplings a^, if the sum in has a 
small number of terms, we will show below that the asymptotic behavior of f+{t) plus the 
D decay data can be used to fix their values, resulting in a fully predictive model for the B 
decay. This procedure assumes that the effect of the heaviest resonances can be absorbed in 
the values of the lighter radial excitations in such a way that the "effective" couplings still 
obey the heavy mass scaling given above. For instance, truncating the sum over resonances 
at 2 = 2, the effect of the heavier resonances can be absorbed, for instance, by 02, defining 



with A2 = M2 — mH and 6i = Mj — M2. As long as the spacing between successive resonances 
becomes smaller, one can neglect 5i/{2mii) corrections. This implies that a|^' has the same 
iriH dependence as 02 so the mn scaling is still valid for the truncated case. The interpretation 
of a|^' is not straightforward. In particular, there is no clear correspondence of this quantity 






(19) 
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with the product of the coupUng g2 times the decay constant of the resonance H2. This, 
however, is not a problem as we will fit the effective parameters to data. Before doing so, 
there is a last theoretical constraint we can impose on f+{t), regarding its behavior at large 
values of 

3.3 Asymptotic Behavior 

The behavior of the form-factor f+{t) for very large values of |t| can be estimated reliably 
in perturbative QCD for exclusive processes (pQCD) In this approach the hadronic 

matrix element is described by the hard scattering transition amplitude folded into an overlap 
integral between the two hadron state wave-functions. To leading order, the hard scattering 
amplitude is approximated by the one-gluon exchange diagrams. The gluon momentum 
satisfies 



where x and y are the momentum fractions of the non-spectator quarks in the initial and 
final hadron, respectively. The positive sign in the last term in (^) corresponds to the 
s-channel process iu Hn, whereas the negative signs corresponds to the t-channel, e.g. 
iH — > UTT as well as to the decay H — > niu. In order for pQCD to be safely applicable we 
need ^ IGeV^. However, the wave-function of a meson containing a heavy quark peaks 
at a: ~ (1 — e), with 



This implies that in the physical region for the decay H iriu the gluon momentum is 
^ IGeV^, with the exception of a negligible large-Q^ tail of the wave- function. This 
casts a serious shadow over the applicability of the one-gluon exchange approximation in 
computing in the physical region for the semileptonic decay, signaling possible large 

corrections. However, outside the physical region and for large enough values of \t\, the 
condition ^ IGeV^ is satisfied. In these two regions, for t <^ and for t much larger 
than the typical mass of heavy resonances, pQCD should yield a very good approximation to 
the form-factor. The known asymptotic behavior of f+{t) constitutes an important constraint 



Q2 = (1 _ xfml + (1 - vYml ± 2(1 - x)(l - y)P ■ 



(20) 




(21) 
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to be satisfied by any calculation. We rewrite the dispersion relation as 

^^(^^ ^ 7 ^ I r^' Im[fl°-'is)]ds ^ 1 fJ^" Im[U{s)]ds 



rn^j, —t) TT Jsa s — t — it 71 JA2 s — t — ie 

^ 1 Im[U{s)]ds ^22) 
TT 7a'2 s — t — ie 



The second term in ( p2D contains the continuum contribution, whereas the third one accounts 
for the region dominated by resonances, for < s < A'^. The last term in (^2]) can be 
calculated perturbatively provided A' is sufficiently large. Its contribution is small for values 
of t in the physical region. At very large \t\, for instance for t <^ 0, the form-factor can also 
be calculated perturbatively. Thus, the asymptotic behavior of (|2^ ) gives 

-1 f 1 z-^' 1 r^'^ 1 

hit) + - im[fl°-'-{s)]ds + - Im[U{s)]ds + p,{t,A')\ 

t [ IT Jso TT JA2 J 

_^ /PQC°(t) (23) 

where the last term in the brackets in (|23|) is the leading asymptotic contribution from the 
last term in (^). The non-perturbative contributions in (^) are, individually, much larger 

than t X /+^'"°(t). Therefore, because /^^'"^(t) is a reliable approximation to the form-factor 
for t — oo, there must be large cancellations among the non-perturbative contributions. 
This leads to a convergence condition or "sum rule" of the form: 

-f-- r /m[/?°°*-(s)]rfs - 1 r /m[/+(s)]ds ~ (24) 

TT Jso TT JA2 

where the equality corresponds to /+*^'"°(t) = —pi{t,A')/t. This sum rule translates our 
knowledge of the asymptotic behavior of /+ (t) as a constraint on the non-perturbative con- 
tributions, which in turn dominate /+(t) in the physical region. In order to actually imple- 
ment this constraint, we can write the integral over the resonant region as a sum over the 
individual resonances. In the narrow width approximation, we have 

1 /-A" Im[f+{s)]ds ^ g, 

2^ JIT— 7 



vr Ja2 s -t - ie ~{ Ml - t 
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where the a^'s are defined in Section 2. In this way the convergence condition now reads 



7-c-^ai~0 (26) 



i=l 



where we defined 

r-A2 



- / /m[/?°°*'(s)]cis (27) 



C = 

Finite width effects are taken into account by the replacement —>■ aili, with the /j's defined 
by 

Ii = -\-- arctan . 28 

The convergence relation of (^) is a very binding constraint for model building. In the next 
section we will explore a specific model for the resonant contribution in order to understand 
the effect of (^6|) on the behavior of f+{t) in the physical region. 



4 Constrained Dispersive Model of f+{t) 

In order to implement the various theoretical constraints from the previous section into a 
calculation of f+{t) we need to specify the resonant contributions. First, we must truncate 
the sum over resonances in (^). As the radial excitations of the H* become heavier, they are 
less relevant to /+(t). On the one hand, heavier resonances contribute with a smaller value 
of TZi even in the narrow width approximation. Furthermore, as finite widths are considered, 
the contributions of heavier and thus broader excitations are additionally suppressed, as can 
be seen in (^). On the other hand, the couplings of the excitations to the ground state, Qi, 
are constrained to obey an Adler-Weisberger sum rule |T^ 

l>g^ + gl + gl + ... (29) 

where gi is the coupling the the z-th radial excitation to ifvr, and additional terms, e.g. from 
orbitally excited states, are not shown. This implies that one cannot add a large number of 
resonances in the cut with large couplings to the ground state. This, together with the mass 
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and width suppression, shows that the truncation of the sum over resonances is a controlled 
approximation. 

In what follows, we will study a Constrained Dispersive Model (CDM) where only the 
first two terms in the sum over resonances are kept. This is partly motivated by the fact 
that only two 1^ radially excited states are known in the light-quark sector. On the other 
hand, the "minimal" choice of keeping only one term will turn out to be incompatible with 
the data on exclusive charm semileptonic decays and the convergence condition (|2^), as we 
will see below. 

The other necessary ingredient to specify the model is the knowledge of the spectrum of 
radial excitations. These resonances ((2S) and (3S) excitations of D* and B*) have not been 
observed in the D or B systems. We will then rely on potential model calculations of their 
masses P]. These models have been very successful in predicting the masses of orbitally 
excited states, and therefore we are confident that the position of the radial excitations in 
the cut does not introduce a sizeable uncertainty. The resulting spectrum explicitly shows 
that the spacings among the IS, 2S and 3S states are, to leading order, independent of the 
heavy quark mass and therefore constitute a property of the light degrees of freedom. We 
take the spectrum of radial excitations to be [^] 

Mf = 2.7 GeV ; Mf = 3.3 GeV . . 

Mf = 6.1 GeV ; = 6.6 GeV ^ ^ 

where the subindex 1 corresponds to the 2S excitation of the H*, etc. In this model, the 
convergence condition (^) now reads 

7 — c — ai — a2~0 (31) 

This condition, together with eqn. (|12D for 7 and the CPTHH calculation for the continuum 



term c defined in (pTl) , leaves only one free parameter in the model. This parameter can 
be fixed by fitting to the observed vr^e^z/g branching ratio. Given that ai/7 and 

02/7 are independent of the heavy quark mass, this procedure results in a prediction for 
B^ TT~i~^i'. The result not only has the correct scaling of f+{t) with the heavy meson 
mass but also implements the HQS properties of the resonances in the cut, namely the scaling 
of the aj's with ttlh- 
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The form-factor emerging in this picture has the form 



mjj.-t\ M| -t ) {Ml- t){Ml - t) 
(M| -s)/m[/^°'^*-(s)] 



(M| - t) TT 7so s - t 



It 



ds (32) 



where we made use of the narrow width approximation. The cutoff in the integral in ( p^ ) 
is given by A ~ {mn + O.TGeV), which corresponds to the maximum center-of-mass energy 
in ifvr scattering for which CPTHH can be used to compute the phase of and at the 

same time coincides with the beginning of the resonance dominated region in the cut. 

In order to fit (|3^ to the D° — > TT~e^Ve data we need values for and g entering in 
the expression (0) for 7. Currently there is only an upper limit on the decay constant from 
D ^ fiu , fij < 0.310 GeV |]TE[. On the other hand, the CLEO collaboration measurement 
of Ds fiu gives /d, = (0.284 ± 0.030 ± 0.030 ± 0.016) GeV 0. We combine this with 



the predictions from lattice calculations |T5[ for the ratio of decay constants, to obtain the 
value /d = 0.24 GeV, to be used in the fit. The H* Hn coupling gets an upper bound 
from the upper limit on the D* lifetime |16|] plus the D* Dn branching fraction [|13 



It is also possible to derive a lower limit on g from D* ^7j. These two combined 

give 0.3 < g < 0.7. We take g = 0.50 for some of our numerical estimates. Finally, 
the prediction for —>■ TT'i^u will depend on the B meson decay constant fs- Lacking 
experimental information on this quantity, we make use of the scaling with heavy meson 
masses and, from the value for we obtain fs = 0.15 GeV. In any event the values of fa, 
Jb and g are external input parameters and not model parameters. Eventually, they will be 
experimentally determined. 



We are now ready to fit (|3^) to —>■ tt e^z/g and fix the value of ai/7. A recent CLEO 



measurement UTS), combined with the value in [|13| for K e^Vf.-, gives Br^^^^{D^ 



n'e^Ve = (3.90 ± 1.5) x 10~^. We do not make use of the more precisely measured — > 
K~l^v mode due to the presence of potentially large SU (3) breaking effects, which we do not 
take into consideration. Thus, at this point, the largest uncertainty in the model prediction 
comes from the experimental uncertainty in the -K^e^Ve branching ratio. It is expected 
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that this will be known to within 5% in the near future . 

With this choice of external parameters, the form-factor of (|3^ ) is shown in the solid 
line of Fig. 2 for the n^i^u mode. Also shown for comparison, are the first term 

in the dispersion relations (dashed-line), corresponding to the S*-pole with 7 given by 
(p!^) (chiral i?*-pole), as well as the prediction from the BSW model (dotted line) pO[. The 
pion momentum distributions for these cases are shown in Fig. 3. As expected, the chiral 
-B*-pole is a good approximation to the full form-factor up to pion momenta of 0{1) GeV. 
Above this momentum, large deviations from the i?*-pole behavior are observed. On the 
other hand, with the choice of parameters made, the model interpolates between the soft 
momentum region, predicted by CPTHH, and the BSW calculation of the form-factor at 



t = 0. This type of procedure was first suggested in [21 1. After all, the relativistic quark 
model calculation of /+(0) in the BSW model constitutes the most reliable prediction of 
quark models for charmless decays of B mesons. Thus, one alternative within our model 
to the need of knowing the external parameters, is to normalize the prediction in B —>■ Triu 



to agree with the prediction from ||20[ at t = 0. As seen in Fig. 3, choosing the central 
values of all the external parameters provides us with a good extrapolation. On the other 



hand, the ISGW2 model is a modification of the non-relativistic quark model of 
which has a harder pion momentum spectrum than the original calculation in order to give 
a better fit to the pion electromagnetic form-factor Ft^{Q'^). Thus the shape of the spectrum 
in the ISGW2 model is a parametrization of the soft physics governing Ft^. The shape of 
the resulting resembles the one obtained in our calculation and is yet another reason 

to believe that the convergence relation captures the correct physics at large values of t. In 
the future, the precise measurement of g, fo and /b will provide an independent prediction 
whithin our model and will help us understand the relation between this simple dispersive 
approach and quark models. 

We now address the possible sources of corrections to this calculation of the form-factor. 
The widths of the resonances are expected to be rather large. When finite width effects are 
taken into account, the form-factor is that given by keeping only the first two terms in 



15 



the sum over resonances, with TZi^2 given by (^) and the convergence relation 

7 - c - aih - 02/2 ^ 0, (33) 

where /i and I2 are defined by (|28D . The widths Fi and r2 are, in principle, not fixed. 
However, they can be estimated by calculating the partial widths to the ground state as well 
as to orbital excitations. To leading order in CPTHH, these interactions are governed by 

2 2 
£ = ^ g,Tr[8,pi^,] + ^ /,Tr [^4^75] (34) 

i=l i=l 

where 7i is defined in (|^) and represents the ground state doublet, £i represent the two radial 
excitation doublets and have the same form as 7i, and the (O"*', 1^) P-wave doublet is given 
by 

S^^-^^[$^,-Sl] (35) 

with Sq and 5*^ the 0"*" and l"*" fields respectively. The couplings gi, g2, fi and /2 are 
unknown. In order to obtain a conservative estimate of the widths we take the values 
91 = 92 = fi = f2 = 0.50. This choice will probably result in an overestimate of the widths 
given it does not respect the condition ([29|). Finally, we account for the the softening of the 
vertices with momentum by the factor 

where k is a typical hadronic scale 0{1) GeV. This results in 

rf ^ (0.25 - 0.35) GeV ; ~ (0.40 - 0.50) GeV 
rf ~ (0.30 - 0.40) GeV ; Tf ~ (0.40 - 0.50) GeV 



(37) 



The resulting prediction for the pion momentum distribution in Ti^i^u is given by the 

dashed line in Fig. 4, together with the narrow width approximation result. The reason for 
the effect being small lies on the fact that, when the widths are incorporated into the fit of 
the model to the —>■ 7r~£^z/ branching ratio, the heavy mass independent parameter ai/7 
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must be slightly modified to upset the width effect. Thus the same type of cancellation takes 
place in the B Triu case. The procedure to obtain a prediction for the B ^ n transition is 
built in a way that is nearly independent of the widths. The model dependence on the cutoff 
A in the continuum integral is also marginally small. This is partly because of the fitting to 
the charm meson decay data, but mostly because the continuum plays a small role overall. 
Finally, we must address violation to the heavy mass independence of the ratios aj/7. The 
heavy mass dependence in these couplings can be written as 



with independent of ttih- Therefore potentially large l/mn corrections partially cancel in 
0^/7 and the use of the flavor HQS in the branch cut is likely to be a very good approximation. 

In Fig. 5 we show the pion momentum distribution for the Ti^i^u decay. The solid 

line represents the CDM prediction, whereas the dashed line corresponds to the chiral pole 
term and the dotted line is the BSW model. The CDM prediction agrees well with both, the 
first term in the dispersion relations and the BSW model. The maximum pion momentum 
is not large enough to cause a disagreement, which arises at ^ 1 GeV and to which only 
the CDM model is sensitive. 

As stated above, our results depend on the external parameters g, fo, fs and Br^'~*'^, 
the TrE'^v branching ratio. In Table I we present the results for the S° tt'E'^v 

branching ratio in units of \Vub\'^ for a particular choice of these external parameters and for 
two sets of masses for the resonances, as well as for the finite width case. The set of radial 
excitation masses we call Set 1 is the one presented in the text, whereas Set 2 corresponds 
to a shift of +100 MeV in the masses. As it can be seen, the result is rather stable under 
these type of modifications. In Table 2 we show results in the narrow width approximation 




7 




Implying 




(38) 
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for fn = 0.24 GeV and for various values of g and Jb- As expected, the predicted branching 
ratio is rather sensitive to the values of these two external parameters, which hopefully will 
be determined in the near future by experiment and/or lattice calculations. 

Although the normalization of the model prediction for Tc^i^u depends on poorly 

known or still unmeasured quantities, one characteristic feature of the model that emerges 
independently of these is the shape of the momentum distribution. As is seen in Figs. 3 and 4, 
at pion momenta above ~ 1 GeV the CDM model predicts that the distribution flattens out, 
as opposed to the characteristic almost-linear growth of the pure i?*-pole behavior. In 
particular, the CDM prediction for the shape at maximum recoil momentum (g^ = 0) is 
consistent with the shape of the perturbative form-factor in the Brodsky-Lepage formalism. 
Considering that the CDM shape is a result of imposing the perturbative behavior in the 
asymptotic region, e.g. for large —q"^, this coincidence is interesting and deserves to be 
studied further. We expect the shape of the pion momentum distribution to be greatly 



constrained by the future CLEO data ||2J]. This will allow to discriminate between the 
CDM approach and the pure -B*-pole model. From Table 1 it can be seen that the intrinsic 
uncertainty of the model is small. This is due to the fact that each change in a model 
parameter leads to a new value of ai/7 to keep a good fit of Br^^'^. This procedure, built 
in the model, makes the predictions for 5° — > Ti^i^u very stable. However, an extraction 
of \Vub\ from the CLEO data is not possible at this point given the large uncertainties in 
external parameters, namely g, fs and Br^^^. 



5 Conclusions 

We have presented an approach to semileptonic form-factors that bridges the gap between 
low and high pion momentum recoil. The dispersion relation described in Section 2 is 
a suitable framework to identify the various contributions. Furthermore, the imposition 
of model-independent theoretical constraints leads to either complete calculations or an 
important reduction in model dependence in calculating the various pieces entering in these 
transitions. 

First, chiral symmetry tells us that , for ^ 1 GeV, the first term in the dispersion 
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relation is a good approximation to f+(t), with 7 given by (|T2|). The singularity closest to 
this region of the H Triu decay is the leading term entering in in CPTHH. The next 
contribution to (Q) comes from the H^*^7r continuum, which dominates the /m[/+(s)] from 
threshold (the branch point) up to a scale where individual resonances start dominating. 
The phase of /+(s) in this region is computed in CPTHH, which is a valid tool up to this 
energy scale. We conclude that the resulting contribution is not the chief modification of 
the H*-pole behavior. 

On the other hand, the radial excitations dominate the Im[fj^{s)] above the scale A, and 
result in large deviations from the H*-pole behavior at values of larger than ~ 1 GeV. 
Their contributions to the branch cut obey a definite HQS scaling with tuh- This imposes a 
constraint that allows us to relate the B -kIv form-factor to that entering in the D TTiu 
mode for all values of the pion momentum. This represents a considerable improvement over 
the simple rriH scaling of the form-factors, which only allows a prediction of the B mode up 
to pion momenta smaller than 1 GeV. 

Finally, requiring that the form-factor has the correct asymptotic behavior as dictated by 
perturbative QCD, we derive the convergence relation or sum rule of (p^ . This translates 
into an additional constraint to be obeyed by the couplings of the radial excitations, as seen 
in dH). 

The model presented in Section 4 has two active resonances in the branch cut, leaving 
only one free parameter, which is fixed by fitting to the observed Br^^'^. The resulting 
form-factor, as seen for instance in (p^ in the narrow width approximation, presents several 
interesting features. Firstly, it respects the constraints from chiral symmetry, HQS and 
perturbative QCD. Secondly, it is a realization of the intuitive idea that the H*-pole behavior 
must be softened by an effective suppression of the coupling 7. In our model this suppression 
is provided by the radial excitations in (|32|) . This can be readily seen if we write the first 
term in (|32D as 

ml,-t[ Mf-t ) - ml.-t[l+v.pJAj ^""^^ 
where Aj = Mj — uih is the gap between the ground state and the i-th excitation, and 
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provides the scale of suppression which in this case is of about (0.8 — 1) GeV. Therefore the 
resonances in the cut give the scale for the effective suppression of the H*-po\e behavior at 
large pion recoil. The model presents a natural explanation for the softening of 7 p5[] . 
It also provides us with a simple rule for the validity of pole dominance in semileptonic 
transitions in general: nearest pole dominance is a good approximation for values of the 
hadronic recoil energy small compared to the gap between the first excitations of the pole 
and the ground state. This is the reason why pole dominance is a good approximation in 
D semileptonic decays as well as in i? ^ D^*Hu decays, but not in i? — > niu, where large 
portions of the rate come from v.p.,^ > Aj. 

As discussed in Section 4, the intrinsic uncertainty of the model is small. This implies 
that a determination of \Vub\ from the measurement of the B —>■ niu branching ratio will 
only be limited by the experimental precision with which Br{D iriu), g, Jd and fs are 
determined. 

The shape of the pion momentum distribution is a very distinctive feature of the model, 
as can be appreciated in Figs. 3 and 4. The comparison with the i?*-pole as normalized in 
CPTHH confirms that pole dominance is a good approximation for ~ 1 GeV. However, at 
larger values of p-,^ the model deviates largely from the characteristic i?*-pole distribution. For 
the same normalization at low p-j^, the branching ratio is about a factor of two smaller than if 
the pure -B*-pole is assumed. The flat distribution at large pion momentum is a consequence 
of the cancellations that result when imposing perturbative QCD as the asymptotic behavior. 
Interestingly, this shape resembles the one that is obtained by computing the form-factor 
directly in perturbative QCD for exclusive processes, characteristic of a dipole fall-off with t 
||26|| , rather than a monopole. Although, as we argued in Section 3.3, the pQCD formalism 
is not expected to give the right answer, this coincidence in shape is an important point to 
bear in mind in an attempt to fully match the non-perturbative and perturbative regimes. 
Present and future CLEO measurements of the pion momentum distribution will test this 
aspect of the model, independently of the values of external parameters. 

Another important test of this approach and the natural next step, is the prediction of 
B piu. In this case the constraints from chiral symmetry are not so clear. However, 
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the HQS and perturbative QCD are still present and largely define the behavior of the 
form-factors, as in the B — > 7r£z/ case. Understanding the B piu to B ^ niu ratio of 
branching ratios, as well as the polarization in the vector mode will turn into a definite test 
of this approach and a major step towards the extraction of |T4b| from exclusive decays in 
the coming B factory era. 
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Tables 





Br{B^^7i-i+u)/\V^,^ 


NWA (Set 1) 


14.00 


NWA (Set 2) 


14.40 


FW (Set 1) 


13.40 



Table 1: CDM prediction of the Br{B^ n-i+u) for g = 0.50, fo = 0.24 GeV and 
/s = 0.15 GeV. Set 1 corresponds to the spectrum of radial excitations of (30) , whereas Set 2 
is obtained by a common shift of +100 MeV. The two first results are in the narrow width 
approximation. The third one includes the effects of the widths as discussed in Section 4, 
with widths given by (37). 
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fs = 0.15 GeV 


/b = 0.20 GeV 


0.30 


19.50 


35.00 


0.50 


14.00 


25.00 


0.70 


7.30 


14.00 



Table 2: GDM prediction of the Br{B^ — > n'i'^u) in the narrow width approximation, for 
fo — 0.24 GeV and with the Set 1 of radial excitation masses. The values chosen for g cover 
the allowed region. The results are shown for two typical values of fs- 
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Figure Captions 



Figure 1: Contributions to the imaginary part of the H niu form-factor f+{t). The 
squares represent the action of the weak current and the circles imply strong interactions. 
The dashed lines are pions and the solid lines are pseudoscalar heavy mesons H unless 
indicated otherwise. 

Figure 2: The form-factor f+{t) as a function of the momentum transfer t. The solid line is 
the CDM prediction for g = 0.5, fo = 0.24 GeV, /s = 0.15 GeV and 5r^^^ = 3.9 x lO'l 
The dashed line is the i?*-pole prediction as normalized by CPTHH in (|I2|) (chiral pole) 
and corresponds to the first term in the dispersion relation (0). The dotted line is the BSW 



model prediction of Ref. |20 



Figure 3: The pion momentum distribution in units of \Vub\'^ as a function of the pion 
momentum. The caption is the same as in Fig. 2. 

Figure 4: The pion momentum distribution in units of \Vub\'^ as a function of the pion 
momentum. The solid line is the CDM prediction in the narrow width approximation. The 
dashed line is the result taking into account finite widths, as given by (37). 

Figure 5: The pion momentum distribution in 7i^i~^i> as a function of the pion 

momentum. The solid line is the CDM fit to Br^^'^ = 3.9 x 10^^ for g = 0.5 and fn = 
0.24 GeV. Also shown are the chiral pole (dashed line) and the BSW prediction (dotted 
line) . 
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